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O . Abstract 
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We analyze the energy density fluctuations contributed by scalar Gelds $ with vanishing 
expectation values, ($) = 0, which are present in addition to the inEaton held. For 
^ . simplicity we take $ to be non-interacting and minimally coupled to gravity. We use 

normal ordering to dehne the renormalized energy density operator p, and we show 
that any normal ordering gives the same result for correlation functions of p. We hrst 
consider massless fields and derive the energy fluctuations in a single mode k, the two- 
point correlation function of the energy density, the power spectrum, and the variance of 
the smeared energy density, {dj^) . Mass effects are investigated for energy fluctuations in 
single modes. All quantities considered are scale invariant at the second horizon crossing 
(Harrison-Zel'dovich type) for massless and for unstable massive fields. The magnitude 
of the relative fluctuations Sp/ptot is of order (ifi/Mpi)^ in the massless case, where Hj 
is the Hubble constant during inEation. For an unstable field of mass m$ <C Hi with a 
decay rate the magnitude is enhanced by a factor ^/m^JV^. Finally, the prediction 
for the cosmic variance of the average energy density in a sample is given in the massless 
case. 

PACS numbers: 04.62. +v, 98.80.Cq, 98.65.Dx 
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1. Introduction and Conclusions 

In the early universe gravity acting on a quantum matter field [1] amplifies energy 
density fiuctuations above the unavoidable vacuum fluctuations [2,3]. These primordial 
fluctuations may serve as seeds for the subsequent formation of the observable large- 
scale structure [4]. 

In the standard literature [5] the fluctuations in the energy density of the inflaton 
field, which drives inflation through the energy density of its vacuum expectation value, 
were thoroughly analyzed. 

In a recent paper [6] we have started to investigate the energy density fluctuations 
contributed by scalar flelds $ with vanishing expectation values, ($) = 0, which are 
present during inflation in addition to the inflaton fleld. For simplicity we take ^» 
to be non-interacting and minimally coupled to gravity. The emphasis of the flrst 
paper has been on the quantum field theoretical concepts and methods. (They are 
briefiy reviewed in the next section.) The present paper contains the application to 
cosmology. We evolve the energy density fiuctuations through the radiation era up to 
the time of matter and radiation equality. A related investigation for the case of classical 
cosmological perturbations can be found in ref. [7] . 

The energy density p of the (non-infiaton) field $ is bilinear in the fiuctuating 
field, and p is therefore not a Gaussian variable. When $ is non-interacting and in 
external gravity, $ is a Gaussian variable and p is x^-distributed. This is in contrast 
to the Gaussian energy density fiuctuations for the infiaton field, which arise from 
the interference between the background field and the fluctuating part of the inflaton. 
Various observable measures of non-Gaussianity in our model are discussed in ref. [8] . 

We use normal ordering to deflne the renormalized energy density operator, and 
we show that any normal ordering prescription gives the same results for correlation 
functions of p. This is so because the difference of two normal orderings is a c-number, 
which drops out in the correlations. 

The cosmological model is given as an inflationary universe, represented by de 
Sitter space with Hubble constant iJi, followed by a radiation dominated universe. The 
total energy density present is always critical, ptot = Pcrit • The transition between 
the two eras can be approximated as instantaneous, because the physical transition 
time is much shorter than the characteristic time for the evolution of the cosmologically 
relevant modes, which have Aphys > e'^'^H^^. The scalar field $ is assumed to be non- 
dominant, therefore it evolves in the background curved space-time. Back reaction of 
$ on the geometry is neglected, and gauge ambiguities are eliminated. We consider 
massless and massive fields, the latter either stable or decaying into radiation. In this 
paper we have not put in a decoherence rate, i.e. the fields $ stay coherent until they 
eventually decay. With respect to decoherence properties the situation is similar to axion 
models (axions, however, have a non-zero expectation value) [9]. Since we have not yet 
included the transfer of the fiuctuations from the resulting massless fields to matter 
which is dominating today (cold dark matter and baryons), we give our predictions at 
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the end of the radiation era. 

The state of the quantum field $ is the Bunch-Davies state initially, during in- 
flation: Every observationally relevant mode, which has today (Aphys)o < H^^, had 
at early times in inflation R/k-p^ys^ — > and is taken to be initially in the Minkowski 
vacuum state. A physical discussion of this state can be found in ref. [6] . 

In section 3 we analyze the contributions of massless fields. We present several 
quantities which characterize the energy density fluctuations: (1) First we derive the 
energy fluctuations SE^: in a single mode k. As long as the mode is super-horizon, 
i.e. between the flrst and second horizon crossing, 5Ej^ grows linearly with the scale 
factor a. This law for massless flelds is independent of the cosmological era. After 
the second horizon crossing SE^j^ ~ a~-^. (2) Next the equal time two-point correlation 
function ^{tj = {5{x)5{x')) of the energy density contrast 5{x) is investigated. At 
the second horizon crossing is independent of the point separation i.e. it is 
of Harrison-Zel'dovich type. On super-horizon scales, for i 3> 'H~^ where Ti is the 
Hubble parameter, the correlation function falls off as For I <^ 'H~^ we flnd a 

logarithmic increase as (log7i£)^ towards smaller scales, ^{tj has a cusp at £ = 2H~^, 
which gives extra strength at large scales. (3) The variance (5^) of the smeared energy 
density contrast 5r is also scale independent at the second horizon crossing. It grows 
as (logTY-R)^ towards smaller scales for R <^ 7i, and it decreases as R~^ on scales 
R ^ H. (4) The power spectrum i.e. the Fourier transform of is such 

that q^P{q) is scale independent at the second horizon crossing, which is again the 
characteristics of Harrison-Zel'dovich type fluctuations. For Ti the quantity q^P{q) 
is ^-independent, while on super-horizon scales {q <^ Ti.) the power spectrum P{q) itself 
tends to a constant. (5) The correlations between the pressure P and energy density 
p are computed during the inflationary era. On super-horizon scales (where they are 
independent of the cosmological era) they are related as 5P = — |(^P- Using the covariant 
conservation of the energy-momentum tensor this again gives the a^-growth of 
for flxed comoving length. 

The magnitude at the second horizon crossing for the relative energy density fluc- 
tuations dp/p is of order (i^i/Mpi)^. Let us assume that during inflation there were 
massless scalar flelds flelds present and take for simplicity A^ = 100 (in the minimal 
supersymmetric standard model we have N = 44). In order to obtain an amplitude of 
order 10~^, we need Hi/Mpi 10~^, i.e. Hi ^ 10^^ GeV is required. For the vacuum 
energy density during inflation, pi := dHjMp^/STT, this 2 ■ IQi^ GeV. 

This is in conflict with the upper bound of pi^/^ ^ 5.2 ■ 10^*^ GeV which was derived in 
[10] from an estimate of the distortions in the cosmic microwave background radiation 
caused by inflationary gravitational waves. 

Mass effects for m <^ Hi are analyzed in section 4 for the energy fluctuations 5E^ in 

— * 

a single mode k. After the super-horizon mode has become non-relativistic, /Cphys ^ 
we flnd a growth of SEj^ ~ a^, i.e. it is enhanced over the SEj^ ~ a law for massless 
modes. The growth of SEj^ stops when Ti. drops below m, i.e. when the Compton 
wavelength Ac = 27r/m becomes sub-horizon. If the massive fleld is stable on the scale 
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of the Hubble time H^q, SEj: is enhanced over the contributions from massless fields by 
a A;-dependent factor at the second horizon crossing. On cosmologically relevant scales 
the enhancement factor is 10^^ or more. To compensate one must take a lower Hi. For 
an unstable massive field which decays with a rate F into radiation the magnitude of 
5Ei^ is only moderately enhanced by a factor ^JrajV independent of the wave number 
k. 

The topic of section 5 is the cosmic variance of p. In a finite sample S", i.e. a finite 
patch of the universe, the spatial average of the density pg — Es/Vs is one measurable 
number. In a quantum field theoretic model one can predict the quantum fluctuations 
of the operator Es, i.e. the variance {5Esf = (El) - {Es)^, where (...) indicates the 
(quantum theoretical) ensemble average. This variance within the ensemble is identified 
with the variance which would be observed when measuring Es in many different patches 
(of size S) within one universe. For a patch size £s ~ 30 /j."^ Mpc there are many patches 
today, and such a variance is accessible to observations and can be compared to models. 
For is — 200 Mpc is (in principle) one measurable number today, but one might 
be able to measure its variance in a decade or so when observations reach much deeper 
than 200 Mpc. For £s — Ho~^ observations can give one single number for pg, 
but its variance can be inferred from COBE measurements under the assumption of 
adiabatic fluctuations together with a modest extrapolation from £ = 2 multipoles to 
£ = 0. This gives the cosmic variance of Sjyg/p ~ 10~^ for is — Hq^. Our model gives a 
contribution to 6'ps/ptot which drops off as l/^/Vg for large sample volumes Vs. Finally 
we relate the experimentally accessible correlation function for {p{x) — ps), which by 
deflnition must have at least one zero, to the quantum fleld theoretic correlation function 
for {p{x) — (p)), which in our model turns out to be strictly positive. On small separation 
scales, £ is, the difference between the two correlation functions is ^-independent 
and equal to the cosmic variance. 



2. Model and Methods 



The background is given by a Friedmann-Robertson-Walker (FRW) space-time 
with spatially flat sections, i.e. we take the mean energy density in the universe to 
be critical, ptot = Pcrit- Using conformal time r], we have ds'^ = a{r])'^ {drf — dx^). 
The Hubble constant during inflation is denoted by Hi. We fix the coordinates in 
such a way that at the time rji of the transition from the infiationary to the radiation 
dominated stage we have 771 = H^^ and a{r}i) = 1. This fixes the scale factor a{r}) for 
our cosmological model as 

^'^ {2-Hiri) ^ ' (2.1) 

a{r^) = Hin (rj > Hf'). 

The Hubble parameter H — a~^9^a at time rj during the radiation era reads 7i = a^^Hi. 
The cosmological model (2.1) is determined by one single model parameter, the Hubble 
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constant Hi during inflation. 



We investigate the energy fluctuations contributed by a neutral scalar quantum fleld 
^» with minimal coupling to gravity. $ has no self-interactions and zero expectation 
value, ($) = 0. Its action reads 



(2.2) 



Variation of the action 5"$ with respect to the field $ gives the equation of motion 
(Klein-Gordon equation) 



+ 2^ — 



d^ 



dr] 



a drj dx^ 



I 2 2 

+ a m 



$ = 0. 



(2.3) 



The energy density p is always defined to be the one measured by a comoving observer, 
who has the 4-velocity u^^ = (a~^, 0, 0, 0) in conformal coordinates. Thus we obtain 
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P = 



2a2 



(a^$)' + (a5f$)' + aW$' 



(2.4) 



The energy-energy correlation at equal times is defined by 

c{x,x') := {n\p{x)p{x')\n) - {n\p{x)\n){n\p{x')\n). (2.5) 

The quantum state \Q) is initially (during inflation) the Bunch-Davies state: Every 
observationally relevant mode, which has today (Aphys)o < Hq^, had at early times in 
inflation R/kph.ya^ and is taken to be initially in the Minkowski vacuum state. Since 
\Q) is translation invariant, C depends only on the physical distance £ = a\x — x'\ of 
the two points. 

We use normal ordering to deflne the renormalized energy density operator 
p{x) •.='N[p{x)]. The key observation [6] is that different choices for the normal ordering 
prescription give identical results for the energy autocorrelation function C{x, x'). This 
is so because two normal orderings N[p] differ by a c-number, which drops out in the 
connected two-point function C{x,x'). 

The mode functions in a spatially flat FRW space-time are eigenfunctions of the 
comoving wave vector k. Because of the translational invariance of the gravitational 
field in 3-space different fc's decouple, and A; is a conserved quantity. We make the 
ansatz for the modes of fixed k 

cpkiv,x) =: Uk{v)e''^^. (2.6) 

The modes Uk{r]) are normalized solutions of the dynamics throughout time, i.e. they 
are solutions of the Klein-Gordon equation (2.3) in both eras of our cosmological model, 
and Uk and drjUk are continuous at the transition. We call Uk, in those modes which have 
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the time dependence e~*'^^ at early times, as kr] — > — oo in inflation. The Uk^iniv) ^-re 
those evolving modes which approached Minkowski single particle waves at early times 
(when their physical wavelengths were much smaller than the Hubble radius). 

Since the results are independent of the normal ordering, we can make the most 
convenient choice, the one adapted to the quantum state \Q). We use the mode expan- 
sion of the fleld operator $(77, x) as 



^{rj, x) 



J (27r 



r 



ak, inUk, iniv) e*^"^ + Ofe, in Uk, iniv)* « 



* „—ikx 



(2.7) 



The state \Q) of the quantum fleld is annihilated by all the operators ak,in, 

O'k, in\^) = 0, for all k. 



(2.8) 



A key tool of ref. [6] was that the computation of the correlation function C{x,x') 
in a given quantum state \Q) can be reduced to the computation of the Wightman 
function of this state. 



w{x,x') := {nmx)^x')\n) = 



d^k 



"^k, in (^)«fe,in(V)*e 



/\* ik(x—x') 



(27r)3 

In the massless case, the two-point correlation function C{x,x') is given as 



(2.9) 



a, 13=0 



(2.10) 



where primed derivatives act exclusively on the primed coordinates. 

The inflationary era was treated in ref. [6]. In particular, for m = the modes 
Uk in have the functional form 



uu,Uri) = -^{i-kf,)e-^^^^+^^^ {m = 0, r,<H-'), 



(2.11) 



where i] = i] — 2H^ ^, i.e. they have the time dependence of a Hankel function, 
Uk, in ~ ■'7^^^ H3/2*^^-*(^^)- The two-point correlation function C was found as 



C{i) 



(27r)^ 



24 UH? 3Ht 



+ 



+ 



£6 2£4 



(m = 0, r]< i^f^). 



(2.12) 



The energy density contrast contributed by the scalar fleld under consideration is 

1 



5{x) 



Ptot 



(p(^)-(p)), 



(2.13) 



where ptot = 37i^Mpj/87r. The two-point function of d{x), 



m = {n\6{x)6{x')\n), 



(2.14) 



is a measure of the relative energy density fluctuations. The two-point functions (2.5) 
and (2.14) are related as ^(£) = C{t)/pl^^. 



3. The Energy Fluctuations for Massless Fields 



3.1. The evolved modes 

After the transition to the radiation dominated era, at times rj > rji = H^^, the 
Klein-Gordon equation (2.3) reads 



'^^ , ^ ^ I ; 2 , 2 tt2 2 

-r-^ H r + k +m HtT] 

ar]'' T] ar] 



(3.1) 



We choose as fundamental solutions for m = 



2kHir] 



(3.2) 



and its complex conjugate v^. They have the time dependence rj }i^^^'^\kr]). The 
mode functions ttfc,in, which correspond to the initial Bunch-Davies state, are linear 
combinations of the two fundamental solutions Vk and v^, 



Uk, iniv) = akVkiv) + PkVkivY 



(3.3) 



The ak and in eq. (3.3) are called Bogolubov coefficients. They are fixed by requiring 
a continuous transition of the mode functions Uk, iniv)j given by eqs. (2.11) and (3.3), 
and their first derivative at time r] = r]i = H^^ . We obtain 



Q!fe = 1 + 



Hi 



k 2fc2 



and 



= 2P- 



(3.4) 



This satisfies the general condition for Bogolubov coefficients Ictfcp — \Pk\'^ = 1- We get 
the functional form of the modes Uk, in during the radiation era as 



Uk, iniv) 



1 



2k Hm 



(3.5) 



There are two types of modes, those which never went beyond the Hubble radius, i.e. 
\hys < 'H~^ at all times, and those modes which spent some time with Aphys > 'H~^. 
The separation point is k = Hi. From now on we look exclusively at the cosmologically 
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relevant modes, which spent a long time outside the Hubble radius, i.e. at modes with 
k < Hi, 

Uk,Uv) = ^^^z-^ {ioik<^Hi). (3.6) 



V2k k kf] 

Since krj = kphys/H we see that well before the second horizon crossing, fcphys -C H, 
Uk,in{v) is time independent. After the second horizon crossing it oscillates in time. 

3.2. Energy fluctuations in one mode 

We start our investigation of the energy fluctuations by looking at one single mode 
k. The mode decomposition (Fourier transform) will be defined in a comoving cube 
with periodic boundary conditions. We compute the fluctuations of the total measured 
energy E^^ within the periodicity cube, contributed by the given mode k, 



{5E^f := {n\El\n)-{n\E^^\n)\ 
The energy of a single-particle quantum measured at a given time is 



(3.7) 



(3.8) 



and the fluctuations of the number of ^-quanta at a given time are 5n^ = 5E^/e^. 

The operator Ey of the total energy in the physical volume V = (aL)^ relates the 
energies Ej^ to the energy density p, eq. (2.4), 



Ev = d^x p{x) = 



(3.9) 



Using the discrete version of the mode decomposition (2.7) of the field operator <^{x) 
we find the operator E^ in terms of bilinears of the annihilation and creation operators 
o/c, in and Ofe^ in^, accompanied by U]^^ in and u^. -^^ with their derivatives. The fluctuations 

of E^j^ are given by — and we obtain 



{5E^? 



a 
1 



k 



where 



(3.10a) 



(3.10b) 



Wk := Uk,in{v)uk,in{V )* 

is Fourier transform of the Wightman function (2.9). 

During inflation the mode functions Wfe, in are given by eq. (2.11), and we obtain 



4k 



phys 



(3.11) 



In the radiation era the mode functions Wfc, in are given by eq. (3.6), and we obtain the 
energy fluctuations 

For very smaU and very large wave numbers compared to the Hubble parameter this 
simplifies to 



_Hl_ 

4^phys 

r2u2 



(for /Cphys -C H) 

(for ?^ « fcphys « — ). 

Let us summarize the time evolution of 5E^ for a fixed comoving wave number k <^ Hi, 
see Bg. 3.1: (1) Before the first horizon crossing early in the infiationary era, SEj^ is 
constant, and Snj^ — > for /cphys ^ Hi. (2) At the first horizon crossing 5E^ = 0{Hi) 
and Srij^ = 0{1). (3) Between the first and the second horizon crossing, when A;phys <S Hi 
and Ajphys <C 'H, the fiuctuations SE^ grow in time linearly with the scale factor a, 
and Suj^ grows as a^. This law remains valid as long as the mode is super-horizon 
sized, independent of the cosmological era. (4) After the second horizon crossing, for 
fephys 3> Tt during the radiation dominated era, the fluctuations SEj^ decrease in time 
~ a~-^. The fluctuations 5n^ do not evolve in time any more. The total growth factor 
for Sn^ between the first and second horizon crossing is Hf /k^, a factor H^/k^ between 
the first horizon crossing and the end of inflation, another factor H^/k^ during the 
radiation era up to the second horizon crossing. 

At fixed time in the radiation era we find 5E^ ~ k~^ for sub-horizon modes, 
fephys ^ 'H. For super-horizon modes, with /cphys ^ H,, we have ~ k~^ . The two 
asymptotes meet at A;phys = Ti.- 

3.3. Two— point correlation function 

— * 

We now sum up the contributions of all the fc-modes to the energy density fluc- 
tuations. To do this we approximate the mode functions Uk, in{v) again by cq. (3.6) 
alone. Thus we will obtain an approximation for ^(£) which is valid only on length 
scales £ ^ aH^^ . But it will fail for the short-distance behaviour of ^(^), on scales i 
which never went beyond the Hubble radius. However, the short-distance behaviour of 
^(£) is universal, independent of quantum state and external curved space-time. This 
term was discussed in ref. [6], and it can be read off eq. (2.12). 

The integral (2.9) for W{x,x') is logarithmically divergent in the infrared, but the 
divergent term is independent of the coordinates of the two points x and x' and therefore 
drops out when computing ^{t). We insert (3.6) into (2.9) and obtain for the infrared 
finite terms of the Wightman function 

H^ 1 



W{x,x') 



(27r)2 24?7?7'Zk 



Ai+log|i?iA_+| + A3_ log|iyiA+_| 

(3.14) 

Ai_ log|iyiA__| - A^^log|iyiA++| 



9 



We have used the notatioiiZk: := — and A-1--1- := ± 77 ± 77'. 

First we investigate the correlation function ^{£) on sub-horizon scales, 
aH-^ <^i^n-^. We obtain 

a platcau-like structure on scales well inside the Hubble radius. Here ^{£) increases 
logarithmically towards the smaller scales. On super-horizon scales, TU ^ 1, we find 

a power law decay. In the range of intermediate length scales £ ^ we need the 
exact expression for the correlation function 



^^^^ ~ (sttMiJ _ (48y4 + 144y2 + 72) + ( 48y5 + 



1 



27y3 144y 



5 y2\ 2/^ 1 5 17 1 

72~36y^ ^ V 192y6 ~ 144y4 ^ 576y2 ^ 64 



(3.17) 

We have used the dimensionlcss variable y = ^7i£, further Ci := log |(y + l)/(y — 1)| 
and £2 := log |y^/(y^ ~ 1)1- The function ^{£) is strictly positive for all £, see Hg. 3.2. 

The correlation function ^{£) is finite for all £ > 0, but it is not diff'crcntiable at 
the point £ = 2^"^ This can already be seen in the Wightman function W{x,x'), eq. 

(3.14). The term involving A has a pole in its third derivatives. This cusp singularity 

has its origin in the instantaneous transition from inflation to the radiation era. We 
can mimic a finite duration by introducing an exponential high frequency cut-off e~^^ 
into the Bogolubov coefficient Pk of eq. (3.4), where T is the fastest comoving time 
scale involved in the phase transition (see the discussion in ref. [11]). However, since 
in any case the time scale T is at most one or two orders of magnitude away from the 
Hubble time during inflation H^^, the smoothing of the cusp occurs on an extremely 
small scale at the time of matter and radiation equality. The effect is invisible on the 
plot. 

The correlation function ^{£) depends on £ only via T-C£, i.e. it is time independent 
if evaluated at fixed ^/^^, where is the Hubble radius. At the second horizon 
crossing, when 7i£ = fixed = 0{1), the correlation is ^-independent and therefore of 
Harrison-Zel'dovich type. The magnitude is of order (i/i/Mpi)^. For instance, we 
obtain y^(£ = 7^-1) = 0.098{Hi/Mpi)^. 



10 



3.4. Variance of the smeared energy density contrast 
We introduce the spatially smeared energy density contrast 



J 5{x)Wr{x). (3.18) 



For convenience, we take a Gaussian window function Wn of scale i?, 

1 

Wr{x) = —= e 2fl2. (3.19) 

The variance is given by (5^) = / d'^x J d^yWii{x)Wji{y)^{x,y). Since the equal time 
two-point function ^ depends only on the physical separation £ of the two points, one 
integration is trivial. We obtain 

We first evaluate (5^) for sub-horizon smearing scales, HR -C 1. For sub-horizon 
separations, H£ <^ 1, the two-point function ^{£) varies only logarithmically, see eq. 
(3.15). Therefore the integral (3.20) is strongly peaked at £ = \/6R, and we obtain 

- \/£ iiy iwS ^'''"'^ 

For super-horizon smearing scales, HR » 1, the exponential in eq. (3.20) can be 
dropped and 

The variance is plotted in Bg. 3.3. For HR ^ 1 we have an R~^ fall-off, in 
contrast to the super horizon fall-off for ^{i). The reason for this is simple: Because 
^{£) falls off rapidly for £ 3> 7Y~^, the integral (3.20) is strongly peaked at £ ^ H~^ , 
independent of the scale i?, if i? » H~^. In this limit the i?-dependence is entirely due 
to the prefactor in eq. (3.20), which reflects the fact that (5^) involves three non-trivial 
integrations. 

3.5. Power spectrum 

Up to now we have concentrated on the two-point correlation function in a;-space. 
For many purposes it is useful to have the Fourier transform of C{£), 

C{q) = \ [ d^£ C(£) e-'^^/"" (3.23) 
a-^ J 

(or correspondingly the power spectrum ^(g) = C((i)/ pl^^)^ where a~^l is the comoving 
separation. By inserting cq. (2.10) we obtain in terms of Uk = Uk^iniv) 

C{q) = / T^<^ \d^Uk\''\dr,u^^_^^\'' + {qk-kY\uk\''\u^^_^^\^ 

^ ^ (3.24) 
-2{qk-P)Re [ukidrjUk)* u^^_^^{drjU^^_j;^y] |. 
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Note that the Fourier transform C{q) of the two-point correlation function C{£) involves 
convolutions of the Wightman function in Fourier space, Wj^ :— A;) — 

\uk,in{v)\^ corresponding to the squares occurring in eq. (2.10). The mode functions 
Uk, in are given by eq. (3.6). The resulting power spectrum C{q) is plotted in Bg. 3.4. 
In the limit q ^ the correlation function C{q) tends to the integral of C{i) over all 
space, see eq. (3.23). Since C{£) is positive everywhere, C{q) must tend to a positive 
constant for q ^ 0. We obtain 

Ciq) ^ ^ jd'i C{i) = 4^ ^ for ?phys « n. (3.25) 

We write C{q = in terms of single mode contributions. From eqs. (3.24) and (3.10) 
we obtain in the continuum limit L — > oo 

(2^ (3-26) 

where := 5Ej^ are the fluctuations in the mean energy density = V~^Ej^ 
in the physical volume V = (aL)^ contributed by the mode k. C{q = 0) is a sum of 
squares, and it cannot vanish unless J-Eg = for all modes k. In the standard treatment 
of inflaton fluctuations only the interference terms between the fluctuations 5$^ and 
the homogeneous background are kept. The interference terms cannot contribute to the 
sum of squares on the right hand side of eq. (3.26), and therefore C(g)inflaton for 
9^0. 

In the sub-horizon wavelength region, gphys ^ 'H, we find numerically 

^ ^ ^ i^) ' for ^ « 9phys « -. (3.27) 

When evaluated at the second horizon crossing, where gphys/'H = fixed = 0(1), the 
quantity q^i{q) is scale independent. The energy density fiuctuations are of Harrison- 
Zel'dovich type. For instance, we find numerically g^C(?phys — '^) — \ (Hi/M-pi)^. 

The relative energy density fluctuations for a flxed comoving q grow in time as 
C{q) ~ 0,^ before the second horizon crossing. This is in contrast to the energy density 
fluctuations in the inflaton scenario which are also of Harrison-Zel'dovich type {i.e. scale 
independent at the second horizon crossing), but which show a super-horizon growth 
as C(?) inflaton ~ o^- Therefore C(q') inflaton ~ q for q-phys < H. 

3.6. Two— point correlations of the energy— momentum— stress tensor 

We investigate the equal-time connected two-point correlation functions of the 
energy-momentum-stress tensor, 

Caicdi^) = {^\ T,f,{x) T^{x') \^) - {n\ (3.28) 

where the hatted indices refer to the local orthonormal frame of a comoving observer. 
For simplicity, we restrict ourselves to the inflationary de Sitter era. These correlation 
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functions depend on the separation vector i — (x — x ')phys- The computation goes along 
the same hne as outhned in section 2, and the correlation functions can be expressed in 
terms of the Wightman function W{x, x') as 

1 (3.29) 

For a = 6 = c = (i = 0we recover eq. (2.10). The tensor structure of C'-g-^ can be built 
using T]^^ and ea = ^a/^- To simplify we average over the directions e, and therefore 
the indices of C.g-^ must be pairwise equal. The possibilities are pp-, pP-, PP , and 

— * — * -| — * 

S ■ (S-correlations, where P = -^{Tn + + is the isotropic pressure, and S is the 
energy flow density (or momentum density), S"' = T°". 

The computation goes along the same lines as outlined in section 2. For super- 
horizon separations, £ » , the energy density and pressure correlations are all of 
the order Hf£~'^ with relative magnitudes 

Cpp{e) = ^Cpp{£), Cpp{e) = -\cpp{i) foriyi^»i, (3.30) 

with Cpp given in eq. (2.12). The ratios correspond to an effective equation of state 

bP = -^Sp. (3.31) 



Using this effective equation of state and the conservation law Tt^^-i, = we can 
derive the growth of the super-horizon perturbations. In a small comoving volume V 
of the external FRW space the conservation law is equivalent to dU = —PdV and to 
dp = — 3(p + P){da/a). Writing p — {p) = 5p we obtain 

d{6p) = -3{6p + 6P)—. (3.32) 

The equation of state (3.31) yields dp ~ a~^, the super-horizon growth factor which we 
have found earlier using a different method, see eq. (3.16). 

— * 

For the autocorrelation function of the energy flow density S we find on super- 
horizon scales 

Cssii) ^ - ^ H,e » 1. (3.33) 

The correlation function of the flow density decreases much faster with £ than the 
correlations of the densities p and P. 

In ref. [12] correlations of the energy-momentum-stress tensor were computed 
which are not averaged over directions. In this case one must distinguish pressures and 
energy flows longitudinal and transversal to £. 
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4. Mass Effects 



4.1. Stable massive fields 

— * 

In this section we investigate the time evolution of 6E^ for a mode k oi a "hght" 
massive scalar with m <^ Hi. As in the massless case of section 3, we consider only modes 
that spend a considerable amount of time outside the horizon, i.e. with /cphys ^ Hi at 
the end of inflation. As long as m -C /cphys the mode is ultra-relativistic and behaves as 
if it was massless. It becomes non-relativistic (/cphys ~ ^) only after the first horizon 
crossing. 

The mass m introduces a second length scale, the Compton wavelength m"-*^, 
in addition to the Hubble length of the external background. Recall that in 

the massless case everything depends on the relation of the two characteristic scales 
= /Cphys and Ti. The horizon crossings occur when /Cphys ~ If m 7^ the ratio 
e^/?i, where e^^ = ^phys + is dynamically relevant for the horizon crossing and not 
kphys/'H. In addition, the ratio kph^ys/m determines the transition from the relativistic 
to the non-relativistic regime. A typical scenario is plotted in fig. 4.1. 

During inflation the massive modes are given in terms of Hankel functions [6], 
Uk,in{v) = I \/ TTfj^Hi H,^(^)(fc?7) of order z/^ = | - m'^Hf'^, where f] = rj - 2Hf^. Using 
eq. (3.10) we obtain for super-horizon modes 

SE^ ~ -3^ (for < Hi and m < Hi). (4.1) 

^ '^phys 

Here we consider only masses (m/Hi)^ <^ ^q. Since cosmologically relevant modes have 

kphys/Hi = C(e~^°) at the end of inflation, factors {k^i^ys/ Hi)"^^ are irrelevant in 
eq. (4.1). As long as the mode is relativistic, i.e. when ~ Ajphys, the fluctuations 

grow in time linearly with the scale factor a. In the non-relativistic regime the growth 
is accelerated to 5E^ ~ a"^, since the physical energy ~ m is no longer redshifted. 

During the radiation dominated era the Klein-Gordon equation is given by eq. 
(3.1). For a highly non-relativistic mode, with /cphys -C to, we can drop the k- 
dependence. We choose as fundamental solutions 

and its complex conjugate, Vkij])* . The mode functions «fe, in are linear combinations of 
Vk and 1?^, constructed in such a way that Uk,ia and its flrst derivative are continuous 
at the time 77 = H^^ of the transition from inflation to the radiation era. We obtain 
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In the modes Uk/m and in the foUowing expressions we neglect powers of {k / Hi)'^^ / . 
The energy fluctuations (3.10) are 

' ' = '^*H:;Vm(^lW j (4.4) 

We distinguish two asymptotic hmits 

^ { _ (4.5) 



(er ~ m > H). 



As long as the mode is super-horizon, SEj^ goes on growing in time as independent 
of the cosmological era. After the second horizon crossing, when ~ m » the 
fluctuations 6Ej^ do not evolve in time any more, SEj^ is constant in time. The time 
evolution of SEj^ is plotted in fig-. 4.2. 

The most interesting effect of the mass m is that the energy density fluctuations 
SEj^ grow faster in time as soon as the super-horizon mode becomes non-relativistic, 
when /cphys ~ Tn and -C 7Y. As a consequence, the fluctuations are much bigger for a 
massive than for a massless scalar fleld. It can easily be read off figs. 3.1 and 4.2 that 
the ratio at the time when A;phys = Ti. (i.e. at the second horizon crossing for fcphys) is 

(|£W2 = ^ (f„rW=W). (4.6) 

{dEj^)m=0 K 

The ratio depends on the comoving wave number k, which equals fcphys at the end 
of inflation. The magnitude of this ratio for modes with Ajphys, eq = -f^EQ, i-e- with 
Aphys — 160 Mpc today, is 



- ' - - ' (for/cphys,EQ = iyEQ). (4.7) 



(5Sfe)m=o V ^EQ 1 eV V 10-28 eV 

Note that the dependence on Hi has dropped out from this ratio since for fixed /Cphys,EQ 
the physical wave number at the end of inflation varies as y/Hi- The massive fiuctu- 
ations are enhanced over the massless fluctuations by an enormous factor for physi- 
cally reasonable non-zero masses. For instance with m ^ 10"^ eV the enhancement 
factor is larger than 10^^. Since massless modes with horizon crossing at equality 
have {dEj:)jri=o ~ Hf/H^Q at the time of horizon crossing, see eq. (3.13), mas- 
sive modes give an order of magnitude of fluctuations compatible with observations if 
Gey)2 ■ (m/10-28 eV)^/^ = C(l). For stable massive fields with m > 10'^ 
eV this requires Hi < 10^*^ GeV, i.e. the energy density during infiation has to fulfil 
pi < (3 • 10^^ GeV)l 
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4.2. Unstable massive fields 

In particle physics all massive bosons are unstable. They will decay into lighter or 
massless particles. We consider the following simple scenario: The massive scalar field 
decays with a rate F into ultra-relativistic particles, i.e. the scalar field has a mean 
lifetime F"^. In addition to the two horizon crossings and the time when the mode 
becomes non-relativistic, there is now a fourth important point in the history of the 
mode, namely the time of decay, where T Ri H. We have included this point in fig. 4.1. 
The decay does not alter the spatial distribution of the energy density on large scales. 
The decay products immediately decohere and thermalize, and they will behave like 
classical radiation. Therefore after the decay the fluctuations Sp in the energy density 
of the massless, thermalized decay products decrease in time as in external FRW 
space-time. 

Instead of evolving in time a true measure of the energy fluctuations in the energy 
density of the decay products, like for instance the two-point correlation function, we 
will give a quick estimate for their magnitude. We continue the energy fluctuations SEj^ 
from the time of decay by redshifting with Qj , clS if the massive mode had transformed 
into a massless (sub-horizon) mode at F = 7i. This scenario is plotted in Bg. 4.3. In 
this case we obtain the ratio 



The ratio is scale independent at the reference point fcphys = Ti-, and does not change 
afterwards. Therefore we can conclude that the energy fluctuations 5p in the massless, 
thermalized decay products have the same scale dependence as in the case of a massless, 
fully coherent scalar field. But their magnitude is enhanced over the massless case 
by a factor ymF^ in the decay scenario. This simple estimate indicates that with 
an unstable massive scalar field the scale invariant (Harrison-Zel'dovich) fluctuations 
are preserved, and the small magnitude {Hi/M-p\f' of the massless case can easily be 
enhanced by a factor of ten or more. 



5. The Cosmic Variance 

5.1. The variance of pg 

In a flnite sample S with volume Vs and mean linear extension £s '■= v^T^, i.e. a 
flnite patch of the universe, the quantum fluctuations of the operator Es of the total 
energy, Es '■= Jgd^x p{x), play a crucial role. (In this section we choose a = 1 at the 
time for which the analysis is performed.) For an infinite sample S the fiuctuations of 
Es are negligible compared to the expectation value (Es), but in a finite S the variance 
(SEs/Es)^ is relevant. 

It is crucial to distinguish between spatial averages and ensemble averages {i.e. 
quantum mechanical expectation values). The spatial averages of an operator A will be 




(4.8) 
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denoted by a bar, A, while for the ensemble average brackets will be used, {A). We note 
that the spatial average of the operator p{x) over the sample S, i.e. pg, is proportional 
to the operator Es, 

The ensemble average {p{x)) of the operator p(x) must be ^-independent in a cosmo- 
logical context, therefore (p) = (pg). The variance of the total energy in the sample is 
(SEs)'^ = {Eg) — {Es)'^, and the variance of the average energy density in Vs is 

mr = ifs) - (ps)'- (5.2) 

For a spherical sample of radius R the variance (^p^)^ is the same as Ptot(^fl) ^ ^^P 
hat window function, Wji{x) = 1 for \x\ < R and Wfi{x) = otherwise. 

Consider now the (equal time) autocorrelation function C(£), eq. (2.5). The spatial 
(double) integral of C{x,x') = C{£) over the sample yields the variance of Eg, 



I d'x I d'x' C{£) = {El) - {Es)' ^ F| • {Sps)'. 
Js Js 



(5.3) 



In the limit V5 — ^ oo the boundary effects of S become negligible, and we obtain 
Jg (fx Jg <fx' C{e) ^Vs- Js d^i C{e). Thus we can rewrite eq. (5.3) as 



(SEs)^ ~ Vs / d^e C{£) 



if Vs^oo. (5.4) 



This is equivalent to {Spg)'^ ~ V^^ C{q = 0) for Vs 00. Since C{q = 0) in our class of 
models is a positive constant, the absolute fluctuations SEs increase with \/Vs, while the 
relative fluctuations 6'ps decrease with 1 / \/Vs- In the standard treatment of the inflaton 
fluctuations only the contributions of the interference terms between background and 
fluctuating part are kept, and 6'ps decreases as for large samples. Because of these 
different asymptotic scaling laws the fluctuations in our model inevitably dominate over 
the fluctuations in the inflaton scenario on large enough scales. Our model predicts 
/ d^£C{£) = llHjH/4207T. Dividing by the total energy density ptot we obtain 



6ps nin AH^ (£nY 



Ptot V 105 3M2j 



3 

(^-^y ioT£s:>£n-=n-\ (5.5) 



5.2. Tiie finite sample correlation function 

We define the operator 

^Ps{x) := p{x)-ps, (5.6) 
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the difference between p{x) and its spatial average in the sample S. We consider the 
(equal time) autocorrelation function Cs{() of the operator Ap5, 

Cs{l) := (Ap5(^)Aps(f' )>• (5.7) 

In a cosmological context Cs{() is only a function of the physical separation £ of the 
two points. The (double) integral of Cs{t) over the sample vanishes by definition, 
Jgd^x Jg d^x' Cs{i) = 0, therefore the correlation function Cs{£) must change its sign 
at some value i in the sample (or at several values). The correlation function Cs{i), 
which by definition has at least one zero, is the experimentally accessible one, see ref. 
[13]. But C{£), which in our model has turned out to be strictly positive, contains the 
ensemble average (p), which is unmeasurable (since only one sample S is available). 

For small separations, £ <^ £s, the two correlation functions are related as 

Ci£) ~ Csi£) + i5psf = Cs{£) + ^ j dH C{£) for £ « £5, (5.9) 

their difference is ^-independent and equal to the cosmic variance. For any fixed £ the 
two correlation functions are identical in an infinite sample. In our model the zero of 
Cs{£) moves to infinity as V5 — > 00. 



We would like to thank Christophe Massacand and Slava Mukhanov for many 
valuable discussions and comments. 
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Figure Captions 



Fig-. 3.1: Time dependence of the energy fluctuations SEj^ for m = 

Fig. 3.2: The energy-energy correlation function C{£) at equahty time. Two domains: 
On scales £ <^ Aeq = i^EQ ~ 13(ri/i^)~^ Mpc a logarithmic increase towards smaller £. 
For £ Aeq a fall-off as £~^. The cusp in the curve is at ^cusp = 2Aeq. 

Fig. 3.3: The variance := Ptot(^fl) a function of the smearing width R at 

equality time. Two domains: For R -C Aeq = -^^eq ~ lS{Q,h'^)~^ Mpc a logarithmic 
increase towards smaller R. On scales R » Aeq a fall-off as R~^. 

Fig. 3.4: The energy-energy correlation function C{q) at equality time. Two do- 
mains: On scales q > Heq = il3{Qh'^)-^ Mpc)"^ a fall-off as q-^. For super-horizon 
wavelengths, q <^ Heq, C{q) is constant. 

Fig. 4.1: Time development of the physical length scales. Three relevant scales for 
a stable massive mode: the Hubble parameter 7i, the mass m, and the energy Cphys- 
A fourth scale for an unstable field, the decay rate F. Three important points: the 
first horizon crossing when Cphys — /sphys = Hi, the time when the mode becomes non- 
relativistic, fcphys = and the second horizon crossing when fiphys — m = Ti. For 
an unstable particle a fourth important point where F — Ti. The point /cphys = Ti. is 
dynamically irrelevant for massive modes. 

Fig. 4.2: Time evolution of the energy fiuctuations SEj^ for a stable massive field with 
<m<^ Hi. 

Fig. 4.3: Time evolution of the energy fluctuations 5E^ for an unstable massive 
field with < m <^ Hi and lifetime F"^. Comparison between the fluctuations SEj^ 
for unstable massive (solid) and massless scalars (dashed). At /cphys = 7i we see the 
enhancement. 
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